Current partition: Nonequilibrium Green's function Approach 
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We present a solution to the problem of AC cur- 
rent partition in a multi-probe mesoscopic conductor 
within the nonequilibrium Green's function formal- 
ism. This allows the derivation of dynamic conduc- 
tance which is appropriate for nonequilibrium situ- 
ations and which satisfies the current conservation 
and gauge invariance requirements. This formalism 
presents a significant generalization to previous the- 
ory: (i) there is no limit in the frequency, and (ii) 
it allows detailed treatments of interactions in the 
mesoscopic region. The formalism is applied to cal- 
culate dynamic conductance of tunneling structures 
with and without assuming wideband limit. 
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The problem of electric current partition in a meso- 
scopic multi-probe conductor is a fundamental issue of 
quantum transport. In the familiar DC situation the 
issue is well understood from the Landauer-Biittiker 
formulationd. Under time varying AC conditions, this 
problem becomes more complicated due to the presence 
of displacement current which-is induced by the AC fields 
according to electrodynamics! 2 ]. In this case, one needs to 
know how to partition the displacement current in addi- 
tion to the particle current for each probes of the conduc- 
tor. Without this knowledge one cannot obtain correct 
results of dynamic conductance: the electric current will 
not be conserved unless the displacement current is taken 
into account. There is another fundamental requirement 
in the transport which is the gauge invariance of the the- 
ory. This simply means that the physics depends only on 
the voltage difference thus shifting voltages everywhere 
by the same constant amount should not not alter the re- 
sults. In the AC transport the gauge invariance will not 
be satisfied without taking into account of displacement 
current. 

Quantum transport under dynamic conditions is the 
focus of a number of recent experimentstTLl. Under low 
frequency AC fields the system is near equilibrium, the 
problem of current partition can be analyzed using scat- 
tering matrix theory (SMT). Biittikei: _and co-workers 
have developed a theoretical formalismfl-u based on SMT 
to analyze linear AC transport and derived the dynamic 
conductance G a p(u)). Here a subscripts denotes a probe. 
By including the contribution of displacement current 
and its partition among the probes, the derived dynamic 
conductance guarantees current conservation and gauge 
invarianceti Obviously, displacement current should be 
even more important for higher frequencies and for situa- 



tions far from equilibrium. In this Letter we address the 
problem of current partition far from equilibrium, and 
with this information we derive a dynamic conductance 
expression which conserves electric current and satisfies 
the gauge invariance under nonequilibrium conditions. 

Far from equiAib*uim, the Keldysh nonequilibrium 
Green's functiont2rtZl (NEGF) has been widely applied 
and many problems of great interests have been analyzed. 
However, within this formalism contribution of displace- 
ment current has not been included. In other words the 
current conservation and gauge invariant condition have 
not been satisfieclEj. A direct consequence of neglecting 
displacement current is to predict incomplete results for 
dynamic conductance, such as that for a parallel plate 
capacitor. The fact of violating current conservation 
within the nonequilibrium thepriea has been recognized 
in a number of publicationsEJ'EJlEil The problem is not 
solved because of a lack of knowledge for partition the 
displacement current. 

Let's consider a quantum coherent multi-probe con- 
ductor with the Hamiltonian 

H = y^ y e ka (t)cl a c ka + H cen {d n ,S n } 

hex. 

+ [ T ka,nC\ a d n + C.C.] (1) 

ka,n 

where tka.it) = £/L + qV a coswt. The first term is the 
Hamiltonian of probes where the AC signal is applied; the 
second term is the general Hamiltonian for the scattering 
region which is a polynomial in {dl_rf„} that commutes 
with the electron number operatoni3 N = X^n^n^"' the 
coupling between the probes and the scattering region is 
given by the last term with coupling matrix Tk a ,n- <4 Q 
(cfe a ) is the creation (annihilation) operator inside the 
a-probe, d) n (d n ) is for the scattering region. 

The particle current (conduction current) inside probe 
a can be calculated using the equation of motion for the 
number of carriers inside that probe. In terms of the 
Green's function and-the self-energy, it can be written in 
the familiar formOo in frequency representation (h — 
1), 

— ±Tr [G r (E + u,,E 1 ) ^<(E 1: E) 

+G < {E + uj,E 1 ) K(Ei,E) - V<(E + u,E x )G a {E u E) 
-Z r a (E + u J ,E 1 )G < (E 1 ,E)] (2) 

where the Green's, function and self-energy are defined in 
the usual manneiOo. The double-time Fourier trans- 
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form is defined as 

F{E 1 ,E 2 ) = J dt x j dt 2 F{t u t 2 )exp[i{E 1 tx-E 2 t2)) 

We will he-, interested in the current component linear 
in voltageEJ, hence we simplify Eq.(|^) by keeping the 
first order in voltage. The dynamic conductance coef- 
ficient G c a p{tu) due to the particle current is defined as 

7 aH = E/3 G «/;H^M- G a/3M is J ust thc frequency 
dependent admittance matrix due to particle current. 
Without introducing confusion, throughout the follow- 
ing we shall use notation G with subscripts to denote 
conductance, while that without subscripts (unless oth- 
erwise stated) to denote Green's function. 

In the presence of AC fields, the particle current alone 
is not conserved due to charge accumulation in the scat- 
tering region Q(t), 2-e,— J^/^ ^ 0. It is the total cur- 
rent which is conservecO^^ Q I^(t)+dQ(t)/dt = where 
I d = dQ/dt is the displacement current. In Fourier space, 



(3) 



where 



QH = -$>/ ^ Tr[g<(E + w,E)] V (4) 



is the AC charge accumulation in the scattering region 
and Og is the small-signal component of the Green's func- 
tion G < defined as J^flgV/j = G < — Gf q in the lin- 
ear regimetil, here Gf q is for equilibrium. This charge 
is related to the displacement field D, of the internal 
Coulomb potential via Poisson equationS (more generally 
the potential satisfies Helmhotz equation): V • D(r, u>) = 
Airp(r, lo) with Q = J pdr. Hence the displacement cur- 
rent — i(u}/Aii) J T)dr is just the time derivative of the 
pile- up charge: J2 a ^a = dQ/dt. Thus the total current 
in probe a is P" 1 = + Current conservation means 
V I tot = 

From Eq.(ph, only the total displacement current 
I d (oj) = —iujQ(uj) is known. We thus need to partition 
it into each probe a. To give an example, the easiest sit- 
uation would be the wideband limit where the coupling 
constants T a are independent of energy. Hence, for an 
amount of charge Q in the scattering region, the proba- 
bility of leaking out through probe a is simply T a /T with 
r the total coupling between the probes to the scattering 
region. This gives the displacement current partition for 
each probe: I d = {T Q /T)I d . 

More generally, we shall use two fundamental require- 
ments of the theory to partition the displacement cur- 
rent: current conservation and gauge invariance. Since 
current in probe a is J Q = ^pGapVp, current con- 
servation means ^2 a Gap = 0; while gauge invariance 
means YIb G a p = 0. We partition the total displace- 
ment current I into the contributions from individual 



probes a in the following form, I a = I c a + A a I d , or 
G af} =G c af3 + A a G d where 

— Trlg^E + ^E)} (5) 

is obtained from Eq. The partition coefficients A a 
must satisfy ^ Q A a — 1 to conserve the total current. To 
determine A a , we apply the gauge invariance condition 
and obtain A a = G^ 7 )/(£ 7 G d ). This gives the 

following form for the dynamic conductance 
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Hence by calculating the Green's functions and using Eq. 
(Q) and (||), we obtain the dynamics conductance which 
now satisfies the current conservation and gauge invari- 
ance conditions. We further point out that result (^J) is 
quite general: it is suitable for AC transport coefficient in 
general terms of u>. It can also be applied to the strongly 
interacting system such as the Anderson model in the 
large U limit by including the contribution of displace- 
ment current. For that model the Green's functions G < 
and G r have been obtained by Ng in Ref. 11. 

Result (^) is in agreement with the formula derived 
using SMT, i.e. Eq. (16) of Ref. |[ However the quanti- 
ties involved in Eq. ra) are calculated within NEGF and 
we now derive these necessary quantities. In the mean 
field approximation it is straightforward to obtain, from 
Eqs.(|) and the expression of /q(w), 
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^0a) + 9 T (3 ^0a 

<-° r a G<)S a p] (7) 



where we have used the abbreviation g = g(E) and 
g = g(E + lu) etc. In the above equation, Go and 
£<,r,a are ^ ne e q U iiik r j ulrL Green's function and self- 
energy function; a < ' r,a is the self-energy function due 
to the AC field: a< = (iq/oj) [T a f - T a f ] and o£° = 

(q/u>) [S q — ^o'q ] i where / is the equilibrium Fermi dis- 
tribution function (at time t = — oo) of the leads far 
away from the scattering region. We have assumed that 
there is no DC bias at time t = — oo. The lesser, re- 
tarded and advanced Green's functions g < ' r < a are de- 
fined as = G ' Q G r ' a and g< = G r a< Gg + 
Gq a r a Gg + Gg ct° Gq . We emphasize that while equi- 
librium quantities such as the Fermi function do appear 
in this formalism, they only represent the "initial" con- 
dition: the system was in equilibrium at t = — oo. For 
t > — oo, the nonequilibrium Green's function approach 
describes the full nonequilibrium dynamic process. With 
the above definitions and quantities, the dynamic con- 
ductance G a /s(uj) can be evaluated for a variety of sys- 
tems near or far from equilibrium. 

The wideband limit. As a first example, it is useful to 
derive an analytic expression of admittance in the wide- 
band limit. Let's consider the case of resonant tunneling 
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through a quantum dot with a single energy level Eq . In 
the wideband limit, the steady state Green's function Gq 
is 
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At zero temperature and through some straightforward 
algebra we obtain from Eq. (^) , 



(8) 



where X is the AC conductance of a noninteracting 
(without Coulomb interaction) symmetrical system ob- 
tained in Refs. ^1 22. The noninteracting result does not 
satisfy current conservation and gauge invariance, while 
the formula (||) does: this is due to the prefactor of X 
in Eq. (||). One can confirm, at the linear order of AC 
frequency, that the wideband result (Q) agrees with the 
scattering matrix theoryQ if we assume a Breit-Wigner 
form of the scattering matrix. 

Beyond Wideband Limit. As a main result of this 
work, we now present the dynamic conductance beyond 
the wideband limit. For this case the situation is quite 
different and we can only obtain results numerically. 
From Eq.(|) and the expression of given above, we 
also note that due to the energy dependence of T, in 
general the dynamic conductance depends on the Fermi 
function in more complicated fashion than simply df /dE. 
In the following let's calculate the admittance of a sym- 
metric double-barrier resonant tunneling structure which 
consists of two (5-function barriers with the same strength 
V . We solve Eq.(H) numerically using a method devel- 
oped by McLennan et. aB$. 

Fig. (1) plots G a g(E,w) as a function of the elec- 
tron Fermi energy E for three values of u. Both real 
and imaginary part of G21 show three peaks due to reso- 
nance transmission through the quantum well. For very 
small u>, e.g. u> = 0.011 (solid line), the real part of 
G2i(E,cu) essentially coincides with the DC conductance 
(not shown); and the imaginary part resembles the emit- 
tance obtained from SMTt3. Substantial deviations is 
observed from these known limits when w is larger. A 
finite u> tends to smear out the resonance behavior of the 
real part of Gi\. In Fig. (2) we plot G2i(Ef,uj) versus 
lj for Ef ~ 0.15 which is near the first resonance en- 
ergy of Fig. (1). The extra peaks of Fig. (2) are due 
to photon assistant tunneling: the electron can absorb 
a photon of appropriate energy and exit from the tun- 
nel structure at another resonance level. This picture is 
confirmed by the inset of Fig. (2) where the energy dif- 
ference and the frequency difference between resonances 
is found to be roughly linear functions of each other, with 
some deviation from perfect linear behavior coming from 
uncertainties due to the width of the peaks. Such a lin- 
ear behavior has been observed in the experiment of Rcf. 
[?]. Finally, in Fig. (3) we plot G2i{Ef,u) versus uj for 
double-barriers with infinite height: for this case there is 



no DC current which can flow through and the system 
becomes a parallel plate capacitor. Due to displacement 
field, our formalism predicts a non-zero G2i(Ef,uj). We 
note that not only the imaginary part of G21 is non-zero, 
even the real part is non-zero. The latter effect gives rise 
to the charge relaxation resistance and is due to .contacts 
implicitly assumed in the transport formalismscj. The 
imaginary part of G21 for this system, dividing 10, is just 
the w-dependent charge response coefficient. It gives the 
electrochemical capacitance of the system when u> —> 0. 
The oscillatory behavior can be traced back to the retar- 
dation effect of Maxwell equations, namely, the Hehnhotz 
equation describing the characteristic potentiaB'cH has a 
term proportional to to 2 which modifies the value and ef- 
fective sign of the induced charge at large u>, leading to 
an oscillatory charge response. 

In summary, we have provided a solution to the prob- 
lem of displacement current partition for multi-probe 
mesoscopic conductors within NEGF. This allows us to 
obtain dynamic conductance in general terms of the AC 
frequency and the result satisfies the current conserva- 
tion and gauge invariance conditions. While for special 
cases (e.g. the wide band limit) analytic results can be 
obtained, in general the dynamic conductance formula 
can be evaluated with well established numerical meth- 
ods. While the nonequilibrium Green's function formal- 
ism can be applied to situations far from equilibrium, 
near equilibrium it must reduce to the same results as 
obtained from the scattering matrix theory. For the ex- 
ternal response, namely the conductance due to particle 
current, Ref. [l4| has provided this connection. Includ- 
ing the displacement current we have shown, using an 
explicit example of double barrier tunneling diode, that 
nonequilibrium Green's function formula indeed recovers 
the scattering matrix result in the wideband limit where 
the dynamic conductance become a Fermi surface quan- 
tity. 

After the submission of our paper, we became aware 
of a recent pape£] which treats photon assisted tunnel- 
ing in the SMT formulation. In this paper, the internal 
potential has been considered explicitly and the gauge 
invariance is satisfied. 
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Units of G21 is 2e 2 /h, of energy E is eV. We have 
set h = l. 

(2) G2i(Ef,oj) as a function of u>, with Ej — 0.1533. 
Solid: read part; dotted: imaginary part. The 
peaks are due to photon assisted tunneling. Inset: 
a roughly linear relationship between the photon 
energy and the resonance energy. Other system 
parameters and units are the same as those of Fig. 
(!)■ 

(3) G2i(Ef,oj) of a parallel plate capacitor as a func- 
tion of w with E F = 0.1533. The solid (dotted) 
line is the real (imaginary) part. Other system pa- 
rameters and units are the same as those of Fig. 
(!)■ 



FIGURE CAPTIONS 



Fig. (1) The real part of G2i(E,u>) as a function of E for 
three values of u. Solid line: u> — 0.01135; dotted: 
u) = 0.1022; dashed: u> = 0.193. Inset: the cor- 
responding imaginary part. The barriers' strength 
V = 8.81, and the quantum well width is IOOA. 
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Figure (2) 



